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PAPER 


Exploring PtS 04 and PdS 04 phases: an evolutionary algorithm based 
investigation 

Horn Sharina," Vinit Sharma,^ and Tran Doan Huan* ^ ^ 


Metal sulfate formation is one of the major challenges to the emissions aftertreatment catalysts. Unlike the incredibly sulfation 
prone nature of Pd to form PdS04, no experimental evidence exits for the PtS04 formation. Given the mystery of nonexistence 
of the PtS04, we explore the PtS04 using a combined approach of evolutionary algorithm based search technique and quantum 
mechanical computations. Experimentally known PdS04 is considered for the comparison and validation of our results. We 
predict many possible low-energy phases of the PtS04 and PdS04 at OK, which are further investigated under wide range of 
temperature-pressure conditions. An entirely new low-energy (tetragonal PA2/m) structure of the PtS04 and PdS04 is predicted, 
which appears to be the most stable phase of the PtS04 and a competing phase of the experimentally known monoclinic C2/C 
phase of PdS04. Phase stability at finite temperature is further examined and verified by Gibbs free energy calculations of sulfates 
towards their possible decomposition products. Finally, temperature-pressure phase diagrams are computationally established for 
both PtS04 and PdS04. 


1 Introduction 

Sulfation (i.e. metal sulfate formation) of noble metal based 
catalysts has been a serious problem to automotive emissions 
aftertreatment systems.^JEllt is well established that Pd is ex¬ 
tremely susceptible towards sulfation (i.e., the PdS04 forma¬ 
tion) in the highly oxidizing and sulfating environment typi¬ 
cally experienced by the aftertreatment catalysts. Unlike the 
easily formed sulfate PdS04 under catalytically relevant con¬ 
ditions, no experimental evidence is available for the existence 
of PtS04 under any circumstances.^^ Despite being a member 
of the same group of the Periodic Table, an intriguing fact of 
non-existence of PtS04 remains as a puzzle and an unexplored 
territory. A question arises why PtS04 does not exist and what 
makes PtS04 different from PdS04? Answers to these ques¬ 
tions may reveal the underlying reason behind the sulfation re¬ 
sistant phenomena of Pt and, in turn, provide some guidance 
for future design of sulfur resistant catalysts materials. 

An experimental investigation based reaction pathway anal¬ 
ysis suggested that the PdS04 formation is primarily due to the 
interaction between SO3 and metal oxide (i.e., PdO) in the cat¬ 
alytically relevant temperature and pressure conditions.l^Nev- 
ertheless, no PtS04 form ation has been observed under similar 
experimental conditions. A recent first-principles compu¬ 
tation based study suggested that the structure of PtS04 should 
be similar to that of PdS04 while assuming a similar nature 
of metal oxides (i.e. PdO and PtO) of Pd and Pt.ElUsing first- 
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principles thermodynamics we have recently predicted that the 
PdS04 formation is indeed favored even at lower temperature 
pressure conditions; however, the PtS04 formation may be fa¬ 
vorable only at elevated pressure conditions. El This outcome 
points out a direction for further investigations of the Pt and 
Pd sulfates under a wide range of temperature and pressure 
regimes, for which comprehensive information on the possi¬ 
ble structural phases is required. Furthermore, PdS04 is sta¬ 
ble towards decomposition to metal oxide (PdO) and sulfur 
oxides (SO2/SO3) below ^ 650 °which suggests that once 
the stable sulfate is formed, it is difficult to desulfate the cat¬ 
alysts. Unfortunately, such information is missing for PtS04 
and needs an attention. 

Our work is premised on the aforementioned mystery of 
contrasting behavior of Pt and Pd metals towards sulfation. 
We extensively explore the possible low-energy structures of 
the yet-to-be synthesized PtS04 and the known PdS04 using 
evolutionary algorithm-based method Unive rsal S tructure Pre¬ 
dictor; Evolutionary Xtalloraphy (USPEX).li 21 Ml xhe thermo¬ 
dynamic stability of the predicted low-energy structures are 
assessed by the evaluation of Gibbs free energy over a wide 
temperature-pressure range, fully considering the vibrational 
contributions calculated within the harmonic approximation. 
Furthermore, we investigate the stability of the predicted 
structures towards decomposition to their possible products. 
In this work, most notably we predict a tetragonal PA2/m 
structure (no. 84 ) to be the lowest in energy for both PtS04 
and PdS04. Interestingly, we find that the experimentally- 
known monoclinic C 2 /c phase (no. 15 ) of PdS04 is ener¬ 
getically competing with the newly identified PA2/m phase. 
From free energies calculations, we propose the temperature- 
pressure phase diagrams for PdS04 and PtS04, predicting the 
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Fig. 1 Relative energetics and structures of the selected low energy phases of PdS 04 and PtS 04 predicted using USPEX method. Pt (or Pd), 
S, and O atoms are represented by silver, yellow, and red colors, respectively. 


Stable phases of the sulfates at high temperatures and/or high 
pressures. 

2 Methods 

Structural phases of compressed matters can now be effec¬ 
tively predicted and discussed at elevated pressures by many 
state-of-the-art computational methods, mostly at the level of 
first principles.El In this work, possible stable structures of 
PdS04 and PtS04 were searched usin g the e volutionary search 
technique embodied in USPEX code.l ^^ l ^^ l This code/method, 
designed to predict the crystal packing from only a knowledge 
of chemical species, compositions, or the molecular geome¬ 
tries, has met tremendous success in correctly identifying and 
predicting the cryst al structures of various classes of systems 
(bulk crystals,f^ 3 QSl nanoclusters,G 21 2 D crystals, 0 ^ surfaces, ESI 
and recently for polymersESEl). In this work, we explored of 
the low-energy configurational spaces of up to four formula 
units of PdS04 and PtS04 per primitive cell, ie., Z < 4 . Struc¬ 
tures with Z > 4 are not considered, and hence, sets a limita¬ 
tion of this work. 

Our first-principles calculations were performed within the 
framework of density function al the ory (DFT) using the pro¬ 
jector augmented wave method^SEIl as implemented in Vienna 
Ab initio Simulation Package (VASPj.l^^l^ While the general¬ 
ized gradient approximation Perdew-Burke-Ernzerhof (PBE) 
exchange-correlation (XC) functional was used throughout 
this work, the energy ordering of the identified structures were 
confirmed to be invariant with the PBEsol^S and the local 
density approximation (EDA) XC functionals. A basis set of 
plane waves with kinetic energy up to 600 eV was used to rep¬ 
resent the Kohn-Sham orbitals while the Brillouin zones were 
sampled by well-converged Monkhorst-Pack k-point meshes, 
i.e., no less than 7 x 7 x 7 . Convergence in optimizing the 
structures was assumed when the Hellman-Feynman forces 


become less than 0.01 eV/A. 

We calculated the densities of states (DOS) of the identified 
structures by the linear tetrahedron method with Blochl cor¬ 
rections. For examining their dynamical stability, the phonon 
frequency spectra calculated using the finite-di splace ment ap¬ 
proach as implemented in the PHONOPY code.l27l2ll To estab¬ 
lish the stability of the predicted phases at finite temperatures 
and pressures, relevant thermodynamic properties were eval¬ 
uated within the harmonic approximation from the computed 
phonon band spectra. FullProf suite^^was used to simulate 
the X-ray diffraction patterns. 

3 Results and discussions 

3.1 Low-energy structures of PdS 04 and PtS 04 

Our evolutionary algorithm based search for low-energy struc¬ 
tures of PdS04 and PtS04, performed at zero pressure (P — 0 
GPa), returned numerous possible candidates. Eight of them 
(six Z = 2 and two Z = 4 structures), which are lowest in en¬ 
ergy for both PtS04 and PdS04, and their energetic informa¬ 
tion are shown in Fig. Of the two common thermodynami¬ 
cally most stable structures of these sulfates, one is described 
by the tetragonal P 42 /m space group (no. 84 ) while the other 
belongs to the monoclinic C 2 /c space group (no. 15 ). It is 
worth noting that the P 42 /m structure can also be obtained by 
substituting Pd/Pt into the Ag sites of the PI (Z = 2 ) structure 
of AgS04.l^ On the other hand, the C 2 /c structure, which 
was experimentally known for PdS04,EIIis similar to that dis¬ 
cussed earlier by Derzsi et a/. El We then found that the C 2 /c 
structure of PdS04 is higher in energy than the PA-2/m struc¬ 
ture by ~ 8 meV/atom, falling within the uncertainty of DFT 
in calculating energies, while for PtS04, this energy differ¬ 
ence is considerably larger, being roughly 20 meV/atom. The 
orthorhombic Ibam (no. 72 ) and the tetragonal P^/n (no. 85 ) 
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Fig. 2 Simulated XRD patterns of two low-energy structures 
predicted for (a) PdS04 and (b) PtS04 at ambient pressure. For 
PdS04, the XRD pattern of the monoclinic C2/c is shown using the 
available experimental data/parameters .ISl The XRD patterns were 
simulated at Cu Kcc with A = 1.54056A. 
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Fig. 3 (Color online) Phonon band structures (left) and density of 
phonon states g{co) (right, given in arbitrary units), of the PA 2 /m 
phase of PtS04 (top panel) and the C2/c phase of PdS04 (bottom 
panel) at P = 0 GPa. For convenience, bands with imaginary 
frequencies, if any, are shown as those with negative frequencies. 
The dotted lines indicate a numerical error of ~ 0.3 THz typically 
resulted from the translational symmetry breaking while calculating 
the XC energies in the real space. 


structures of both PdS04 and PtS04 are relatively similar in 
energy, residing at > 50 meV/atom above the PA^Itn. Two 
Z — A triclinic structures examined, namely a-Pl and j 3 -Pl 
(both no. 1 ), are about 30 — 50 meV/atom above the PA^jm 
structure. The last two structures, i.e., 7222 (orthorhombic, no. 
23 ) and 74 (tetragonal, no. 82 ), are about 75 — 100 meV/atom 
higher than the PA^jm structure. We note that these structures 
are slightly below the C 2 /c (Z = 16 ) structure recently deter- 
minedl^for AgS04. This energy ordering remains essentially 
unchanged when PBEsol and LDA were used (also see Fig. 
[T]). Crystallographic information of the predicted structures is 
given in the Supporting Information (see Table SI). 

The predicted low-energy structures of both sulfates consist 
of tetrahedral SO4 groups, where O atoms are associated to 
four different SO4 tetrahedra coordinating the Pd/Pt atoms on 
a plane. The local chemistry at the anionic site, the topology 
and the connectivity of the crystal networks are also qualita¬ 
tively similar in both sulfates. In general, these eight structures 
can be classihed into two groups. The first structure type, with 
a non-layered 3 -D network, contains oxygen atoms from the 
SO4 unit which act as a bridge (for example: PA^jm, Cljc, 
74 , and j 3 -Pl phases) linking metal atoms. The second struc¬ 


ture type involves some two-dimensional motifs with isolated 
layers of Pd/Pt and SO4 tetrahedra. The remaining four struc¬ 
tures, i.e., Ibam, P^jn, 7222 , and a-Pl phases, belong to this 
class. 

We further analyzed the selected low-energy structures by 
simulating the X-ray diffraction (XRD) patterns. In Fig. 
(top two panels), we show the XRD patterns simulated for the 
C 2 /c and P42/m structures of PdS04 along with the available 
experimental XRD data of the Cljc phase.® In the bottom 
two panels of Fig. we show the XRD patterns simulated for 
our predicted C 2 /c and P42/m phases. Overall, the simulated 
XRD patterns are in good agreement with the available exper¬ 
imental data. ®The additional simulated XRD patterns (of the 
other predicted phases) are given in the Supporting Informa¬ 
tion S 2 . 

3.2 Dynamical and thermodynamic stabilities 

Next, we examined the dynamical stability of the predicted 
structures of PtS04 and PdS04 using the calculated phonon 
band structures. No imaginary modes exist throughout the 
Brillouin zones of these structures, demonstrating that they 
are dynamically stable. For illustration, we show in Fig. 
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Fig. 4 Gibbs free energies with entropic contributions, G calculated at T = 0 K, T = 300 K,T = 700 K, and T = 1000 K for the identified 
low-energy structures of PdS04 (panel a) and PtS04 (panel b) are shown as functions of pressure P. Data is given by symbols while curves 
are guides to the eye. 


the phonon spectra and the phonon density of states g{(o) we 
calculated for the lowest-energy structures of each compound, 
i.e., the P 42 /m and C 2 /c structures. Similar information for 
all other predicted structures can be found in the Supporting 
Information S 3 . 

The phonon spectra of these structures, calculated at 0 K, 
allow estimating the vibrational contribution /Vib( 7 ’) to the 
Gibbs free energy G{P,V,T) = Zsdft +Fvib{T) + PV within 
the harmonic approximation via 


FvibiT)^rkBT [ d(0g{(0)\n 
Jo 


2 sinh 



( 1 ) 


where, r is number of degrees of freedom in the unit cell, 
is the Boltzmann’s constant, h is the reduced Planck’s con¬ 
stant, and g{co) is the normalized phonon density of state at 
frequency co. In addition, the enthalpy fsopT + PV was cal¬ 
culated by slowly optimizing the investigated structures under 
gradually increasing pressure, starting from P — 0 GPa. For 
hard crystalline materials, the method provides a reliable way 
to access their thermodynamic stability. For hard crystalline 
materials, this method typically leads to an excellent agree¬ 
ment with experimental data.l^ 

The calculated free energies G{P,V,T) are summarized in 
Fig. El suggesting that the PA-i/m phase of PtS04 is ther¬ 
modynamically stable at low pressures. Within this regime, 
the C 2 /c structure of PdS04 (which is experimentally estab- 
lishedl^ is different from the P42/m structure by no more 
than ±2 meV/atom at low and high temperatures. Therefore, 
these two phases are considered to coexist at low pressures. 
The formation of the C 2 /c phase, which is observed even 


at low temperatures conditions, may be driven by kinetics, 
known under the empirical Ostwald’s steps rules in crystal nu¬ 
de ation. 

Both PdS04 and PtS04 undergo several structural phase 
transitions at elevating pressures. For PdS04, the orthorhom¬ 
bic Ibam phase is stable between ^ 10 and ~ 60 GPa before 
transforming to the triclinic j 5 -Pl phase. The Ibma-to-P-Pl 
phase boundary depend very weakly on temperature. Unlike 
PdS04, the C 2 /c phase of PtS04 is thermodynamically stable 
only at high temperature (> 700 K) and intermediate pressure 
(10 — 30 GPa) conditions while the Ibam phase is stable at 
lower temperatures (< 700 K) and elevated pressure (10 — 60 
GPa) conditions. The transition between the Ibam phase to the 
74 phase occurs at roughly around 60 GPa. For both PdS04 
and PtS04, the 1222 phase is unstable over the whole range of 
pressure examined. 

Using the calculated free energies G{P,V,T) (as shown in 
Fig. | 4 |, we constructed the temperature-pressure phase dia¬ 
grams of both sulfates and show them in Fig. The phase 
diagrams display a map of the stable phases over the range 
of T-P conditions. Most importantly, we observed that for 
PdS04, both the tetragonal PA^jm and the monoclinic C 2 /c 
phases coexist at atmospheric pressures while for PtS04, the 
PA2/m phase is the sole candidate at the same conditions. Fur¬ 
thermore, the Ibam phase dominates the 10 — 60 GPa region 
for both cases, which could be an interest of exploration for 
the high pressure applications. 

While experimental studies suggest that PdS04 decom¬ 
poses above ^ 900 K^ estimation of free energy of reaction 
(AG) allows us to evaluate the thermodynamic stability (reac- 





































Fig. 5 Computed phase diagrams of PdS04 (left panel) and PtS04 
(right panel). Thermodynamically stable phases are shown as 
indicated hy their space group symbols. Shaded area indicates the 
coexisting regime of both the PA^jm and C2/c phases of PdS04. 

tion feasibility) of the compound towards the decomposition 
into possible products. The feasibility of a reaction depends 
on the sign of AG, which is equal to A// — TAS, where AH is 
the change in enthalpy and AS is the change in entropy. The 
AG of the reaction can be expressed as: 

n n 

Gproducts frreactants- (2) 

!=1 !=1 

In this work, we considered the decomposition reaction of 
Pd(or Pt)S04 towards their respective most stable metal ox¬ 
ides and sulfur oxide species [i.e. Pd(or Pt)S04 —f Pd(or 
Pt )0 + SO3]. For example, the computed AG values at 300 K 
were ~ —60 kJ/mol and ~ —40 kJ/mol for PdS04 and PtS04, 
respectively. Furthermore, we evaluated the free energy of 
the decomposition of the sulfates to their respective elemental 
species (i.e. Pd (or Pt) SO4 —Pd(or Pt) H- S H- 2O2). The 
computed AG values were in the range of ~ —500 kJ/mol at 
300 K. Results suggest that PdSOd is stable towards decom¬ 
position to PdO and SO3 below 775 K whereas PtS04 stability 
towards PtO and SO3 remains below 650 K. Similarly, PdSOd 
is stable towards the decomposition to the elemental compo¬ 
nents below 870 K whereas PtS04 is stable below 800 K. Fur¬ 
thermore, we computed the dG for the reaction Pd(or Pt)S04 
—> Pd (or Pt)S + 2O2, which further supports the stability of 
the sulfates in realistic temperature (< 700 K) and pressure 
conditions. In general, our calculations show that PdS04 is 
more stable than PtS04 towards decomposition for a particu¬ 
lar temperature. Our results are in good agreement, given the 
computational error range in energetics, with the available ex¬ 
perimental results of PdS04 decomposition stability. Further¬ 
more, synthesis of PtS04 seems feasible in the future given 
the kinetic barriers are easy enough to cross. The free energy 
(AG) versus temperature (T) plot is provided in Supporting 
Information (S 4 , Figure S 3 ). 


To further confirm whether these sulfates are stable or not 
with respect to the pool of all possi ble product species, a linear 
programming (LP) algorithmEMl has been employed. Here, 
a Pd(or Pt)S04 compound is considered to be stable when AE 
(the DFT energy relative to the best outcome from the LP) is 
negative. The energy difference, AE, can thus be written as 

n 

AE = Pd (or Pt)S04 - min ^ aPi, ( 3 ) 

1=1 

where Pj represents all the possible stable chemical species 
(i.e. for PdS04: Pd, PdO, PdS, SO3, SO2, SO, S, and 
O2; for PtS04: Pt, Pt02, PtO, PtS, SO3, SO2, SO, S, and 
O2). For example, the equation for PdS04 becomes AE — 
PdS04-min(ciPd„, -|-C2Pd„2 0o2+C3Pdfo3S3o3 + 

C5S65O205 +C6S65O0J TcySfo, -f C8O208). Then, the LP prob¬ 
lem is solved with the constraints 

= 1 , ^ Vi = 1 , and = 4 , ( 4 ) 

i i i 

where a,, bi, and o,- represent Pt(or Pd), S, and O content of 
a species, respectively. Above constrains ensure the correct 
stoichiometry of Pd(or Pt)S04 and with 

Ci>l, ( 5 ) 

which warrants that only the references containing Pt(or Pd), 
S, or O are taken into account. With all DFT computed ener¬ 
gies of the species, we obtained all the optimized c, and AE for 
each case. Consistent with our free energy of reaction anal¬ 
ysis, negative AE values (i.e. — 0.87 eV and — 0.70 eV for 
PdS04 and PtS04, respectively) were obtained, which con¬ 
firmed the stability of the sulfates. Interestingly, we obtained 
mono-metallic oxides (PdO and PtO in the case of PdS04 
and PtS04, respectively) and SO3 as possible decomposition 
products,consistent with our reaction free energy analysis, and 
unity (as expected) for all c, values. 

3.3 Electronic structures 

We investigated the electronic structures of all low-energy 
phases of PdS04 and PtS04 by computing the total density 
of states. Overall, our results show no significantly different 
behavior between the phases of both sulfates. For better un¬ 
derstanding the DOS can be divided into three main groups. 
First, the lower valence bands (between —2 eV and —4 eV) 
originate due to mixing of the valence d and p states of Pd(or 
Pt) and the O atoms. In this region, Pd(or Pt) {d) bands are 
found to be highly resonant with the O (p) bands. We also no¬ 
ticed that some pronounced mixing between the segments of 
O (p) bands lying above and below the valence Pd(or Pt) (d) 
bands. Second, in the vicinity of the Fermi level the valence- 
band maximum are dominated by Pd(or Pt) (d) states. Third, 
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the bottom of the conduction band consists 3p states of S and 
O (2p) states. Further detail can be found in the Supporting 
Information S5. 

4 Conclusions 

In summary, we explored the mystery related to nonexis¬ 
tence of PtS 04 using first-principles thermodynamics com¬ 
bined with the evolutionary algorithms based method. Our ap¬ 
proach is validated by also studying the experimentally known 
phases of PdS 04 . Many low-energy structures are predicted 
and analyzed for the stability in a wide range of tempera¬ 
ture and pressure conditions. At low pressures, we identify a 
tetragonal PA 2 /m structure (of the AgS 04 type) which appears 
to be the thermodynamically most stable phase of PtS 04 . In 
case of PdS 04 , this phase is predicted to coexist with the ex¬ 
perimentally known C2 /c phase. These sulfates are also pre¬ 
dicted to undergo several phase transitions at elevated temper¬ 
atures and/or pressures. Based on the computed Gibbs free 
energies, we constructed phase diagrams which provide such 
the reliable information about the phases stability, the phase 
transition, and their boundaries up to 100 GPa and 1000 K. 
The phase diagrams confirmed the existence of experimen¬ 
tally observed monoclinic C2/c phase of PdS 04 at the am¬ 
bient conditions; however, this phase may not be seen in the 
case of PtS 04 in similar conditions. Nonetheless, Ibam phase 
remains one of the promising stable phase for both cases at 
high pressure conditions. Both sulfates were stable towards 
decomposition to their possible products well above the room 
temperature, which also suggests the possibility of PtS 04 syn¬ 
thesis in the future. In general, we provide a detailed infor¬ 
mation on the phases and their stability of PdS 04 and PtS 04 
which can be helpful to understand the sulfating nature of Pd 
and design/scan promising new sulfur resistant materials. 
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* Corresponding author: huan . tran@uconn . edu 


References 

1 H. N. Sharma, S. L. Suib and A. B. Mhadeshwar, in In¬ 
teractions of Sulfur Oxides With Diesel Oxidation Cata¬ 
lysts (DOCs), American Chemical Society, Washington, 
DC, 2013, vol. 1132, pp. 117-155. 

2 A. Russell and W. S. Epling, Catal. Rev., 2011, 53, 337- 
423. 

3 H. C. Yao, H. K. Stepien and H. S. Gandhi, J.Catai, 1981, 
67, 231-236. 

4 C. P. Hubbard, K. Otto, H. S. Gandhi and K. Y. S. Ng, 
Catal. Lett., 1994, 30, 41-51. 

5 C. P. Hubbard, K. Otto, H. S. Gandhi and K. Y. S. Ng, J. 
Catal., 1993, 144, 484^94. 

6 S. Koutsopoulos, S. B. Rasmussen, K. M. Eriksen and 
R. Eehrmann, Appl. Catal. A: Gen., 2006, 306, 142-148. 

7 M. Derzsi, A. Hermann, R. Hoffmann and W. Grochala, 
Eur. J. Inorg. Chem., 2013, 2013, 5094—5102. 

8 D. L. Mowery, M. S. Graboski, T. R. Ohno and R. L. Mc¬ 
Cormick, App/. Catal. B: Em., 1999, 21, 157-169. 

9 G. Corro, React. Kinet. Catal. Lett., 2002, 75, 89-106. 

10 U. Kohler and H. W. Wassmuth, Surf ScL, 1983, 126, 
448^54. 

11 H. Sharma, V. Sharma, A. Mhadeshwar and R. Ram- 
prasad, J. Chem. Rhys. Lett., 2015, 1140-1148. 

12 A. K. Neyestanaki, E. Klingstedt, T. Salmi and D. Y. 
Murzin, Fuel, 2004, 83, 395^08. 

13 C. W. Glass, A. R. Oganov and N. Hansen, Comput. Rhys. 
Commun., 2006, 175, 713-720. 

14 A. R. Oganov and C. W. Glass, J. Chem. Rhys., 2006,124, 
244704. 

15 E. Zurek and W. Grochala, Rhys. Chem. Chem. Rhys., 
2015, 17, 2917-2934. 

16 A. R. Oganov, A. O. Lyakhov and M. Valle, Acc. Chem. 
Res., 2011,44, 227-237. 

17 X.-E. Zhou, X. Dong, A. R. Oganov, Q. Zhu, Y. Tian and 
H.-T. Wang, Rhys. Rev. Lett., 2014, 112, 085502. 

18 A. O. Lyakhov, A. R. Oganov, H. T. Stokes and Q. Zhu, 
Comp. Rhys. Comm., 2013, 184, 1172-1182. 

19 Q. Zhu, L. Li, A. R. Oganov and P. B. Allen, Rhys. Rev. B, 

2013, 87, 195317. 

20 V. Sharma, C. C. Wang, R. G. Lorenzini, R. Ma, Q. Zhu, 
D. W. Sinkovits, G. Pilania, A. R. Oganov, S. Kumar, G. A. 
Sotzing, S. A. Boggs and R. Ramprasad, Nat. Comm., 

2014, 5, 4845. 

21 Q. Zhu, V. Sharma, A. R. Oganov and R. Ramprasad, J. 
Chem. Rhys., 2014, 141, 154102. 

22 W. Kohn and L. J. Sham, Rhys. Rev., 1965, 140, A1133- 
A1138. 

23 P. Hohenberg and W. Kohn, Phys. Rev., 1964, 136, B864- 




B871. 

24 G. Kresse and J. Furthmiiller, Phys. Rev. B, 1996, 54, 
11169. 

25 G. Kresse and J. Furthmiiller, Comput. Mater. ScL, 1996, 
6 , 15-50. 

26 J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, 
G. E. Scuseria, L. A. Constantin, X. Zhou and K. Burke, 
Phys. Rev. Lett., 2008,100, 136406. 

27 A. Togo, F. Oba and I. Tanaka, Phys. Rev. B, 2008, 78, 
134106. 

28 K. Parlinski, Z. Q. Li and Y. Kawazoe, Phys. Rev. Lett., 
1997, 78, 4063-4066. 

29 J. Rodriguez-Carvajal, Physica B, 1993, 192, 55. 

30 P. Malinowski, M. Derzsi, Z. Mazej, Z. Jaglii, B. Gawe, 
W. asocha and W. Grochala, Angew. Chem. Int. Ed. (En¬ 
glish), 2010, 49, 1683-1686. 

31 T. Dahmen, P. Rittner, S. Boger-Seidl and R. Gruehn, J. 
Alloy Compd., 1994, 216, 11-19. 

32 M. Derzsi, A. Budzianowski, V. V. Struzhkin, P. J. Ma¬ 
linowski, P. J. Leszczynski, Z. Mazej and W. Grochala, 
Cryst. Eng. Comm., 2013, 15, 192-198. 

33 D. L. Mowery and R. L. McCormick, Appl. Catal. B: Env., 
2001, 34, 287-297. 

34 T. D. Huan, V. Sharma, G. A. Rossetti, Jr. and R. Ram- 
prasad, Phys. Rev. B, 2014, 90, 064111. 

35 I. E. Castelli, D. D. Landis, K. S. Thygesen, S. Dahl, 
1. Chorkendorff, T. F. Jaramillo and K. W. Jacobsen, En¬ 
ergy Environ. Sci., 2012, 5, 9034-9043. 

36 V. Ozolins, E. H. Majzoub and C. Wolverton, J. Am. Chem. 
Soc., 2009,131, 230-237. 






